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Abstract 

We discuss the possibility that high-frequency QPOs in neutron-star binary systems may result from forced resonant 
oscillations of matter in the innermost parts of the accretion disc, excited by gravitational perturbations coming from 
asymmetries of the neutron star or from the companion star. We find that neutron-star asymmetries could, in principle, 
be effective for inducing both radial and vertical oscillations of relevant amplitude while the binary companion might 
possibly produce significant radial oscillations but not vertical ones. Misaligned neutron-star quadrupole moments of a 
size advocated elsewhere for explaining limiting neutron star periods could be large enough also for the present purpose. 
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1. Introduction 

Quasi-periodic oscillations of X-ray brightness (QPOs) have 
been observed in a number of accreting binary systems 
, containing c o mp act objects, both with neutron stars (see 
■ Ivan der KhsL [2000: Barret et al., 200£ j for a revie w) and 
I with black holes (,Remillard fc McClintockl . [200i). They 
, can have low frequencies (Hz) or high frequencies (kHz). 
The observed kHz frequencies are comparable with the 
Keplerian and epi cyclic frequen cies in the inner parts of 
, the accretion disc (jTorokl |2005| ) and the question arises of 
whether they might be associated with forced resonant os- 
cillations of the inner disc material. In order to initiate 
these, some perturbation mechanism would be required. 
Here, we focus on neutron-star systems and investigate the 
possibility that gravitational perturbations caused either 
by the binary companion or by asymmetries of the neutron 
star might provide this mechanism. (The case of the binary 
companion could be relevant also for black hole systems.) 
We note that perturbations coming from the neutron star 
can only be relevant for the innermost parts of the disc 
(because of the rapid fall-off of the force with distance) and 
hence are mainly associated with the picture for kHz QPO 
frequencies rather than with that for the lower frequencies. 
It is necessary that the mechanism should be a resonant 
one since otherwise the response produced would certainly 
be much too small to be of interest. 
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The two types of perturbation (from the neutron star 
and from the binary companion) clearly induce different 
behaviours: the frequency of the varying force arising from 
the influence on the disc of the binary companion is es- 
sentially equal to the disc rotation-frequency, whereas the 
main frequency of the force caused by asymmetries of the 
neutron star is equal to the difference between the rota- 
tion frequencies of the disc and the neutron star (jPetril 
2006). Resonance occurs at those points where one of the 
intrinsic oscillation frequencies of the disc matter coincides 
with the forcing f r equen cy. Following the discussion by 
iLandau &: Lifshit3 ()1976f ) for test particle motion (which 
would need to be modified for fluid elements), the growth 
in amplitude of the oscillations of the particle within the 
linear regime of forced resonance is given by 



where /p is the amplitude of the variations of the force, to 
is the frequency and mp is the mass of the particle. This 
linear regime ends when the oscillation amplitude a(t) be- 
comes large enough so that non-linear phenomena and/or 
dissipative processes become relevant. Note that a{t) grows 
linearly with time in this regime and so can become quite 
large even when the variations in the perturbing force are 
small. 



Z. Stuchlik et al.: Gravitational excitation of QPOs (RN) 



2 




The vertical force oscillates around its mean value with 
amplitude 



Av = (Gm/Rl) cos 9 A 

[1 - 2a; sin 61a + x^) 
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and the radial force oscillates with amplitude 



= {Gm/R 



2\ 2 



(1 — a; sin 9^) (l — 2.t sin 9^ + x 
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These oscillations have an anharmonic character which 



related to it. 



2.2. Binary companion 

The binary companion is taken to be orbiting the neutron 
star at a constant distance d, with angular velocity Q^. 
The vertical and radial components of the perturbing force 
acting on the accreting material are then given by 

Fbv(<) = (GM^/d^) COS 9^ 

X [1 - 2{R/d) sin 9^ cos cu^t + [R/df] (6) 



F^n{t) = (GAfe/rf^) [sin^BCOS^Bt- 

X [1 - 2{R/d) sin 9^ cos oj^t + [R/df] ^^'^ 



(7) 



Figure 1. Schematic picture illustratmg the generation of , . • , i ;, , , ;, ;, 

., ,. , 4. u- f • 4. • 1 4.- means that, when iourier analysed, they show both the 

a gravitational perturbing force m an equatorial accretion , . ^ , , r 

J. , • 1 11 i • )) J.^ c r J.1 i basic frequency uj^ and also some additional frequencies 
disc by a single mountain on the surface of the neutron ^ j a -l 

star (A/a7 Ra) or by the binary companion {M^, R^)- 
2. Gravitational perturbing force 

We consider here the situation illustrated in Fig.[T] with a 
basically isotropic neutron star of mass A/a and radius i?A 
spinning about its rotation axis with angular velocity fl^ 
and with the symmetry plane of the accretion disc being or- 
thogonal to the rotation axis. We use spherical polar coor- 
dinates {R, 9, if) with the origin at the centre of the neutron 
star. Possible neutron-star asymmetry is approximated by 
a point-like source with mass m located on the surface of 
the star at angle 9 — 9j^. The binary companion (with mass 
Mb) is taken to be moving on an orbit at a constant distance 
d from the neutron star at 9 = 9^. For thin discs, the an- 
gular velocity profil e ^d{R) is generally well- approximated 
as being Keplerian ([Novikov &: Thornel . |1973[ ) but for thick 
(or slim) discs, there is a deviati on away from this becaus e 
of the action of pressure forces (jJaroszyhski et ahl . 119801 ). 
We here determine the radial and vertical components of 
the gravitational force produced by the perturbing sources 
in a purely Newtonian way. This is an approximation, but 
we do not expect that a relativistic analysis would greatly 
change the qualitative features of our results. For simplic- 
ity the force will be determined in the equatorial plane, 
i.e., in the symmetry plane of the disc; this is completely 
correct for thin, Keplerian discs, and gives good estimates 
for slim discs. We determine the time evolution of the per- 
turbing force components for a fixed point on the disc with 
coordinates (i?, = 7r/2,(/3 — fJdt), making the restriction 
Ra < i? < 10 i?A since kHz QPOs are being considered. For 
the binary companion, we assume d ^ R, R^. 

2.1. Neutron-star asymmetry 

Here, and in the following, we consider the force acting on a 
comoving unit mass element of the accretion disc at a given 
radius R. Using the quantities x = Ra/R and lja = |^^a — 
fid I , the vertical component of the perturbing gravitational 
force is given by 

FAv{t) = (Gm/Rl) x^ cos 9 J, 

X (l — 2a:: sin ^A coswa^ + 
while the radial component is given by 
^ar(^) — {Gm/R\) x^{l — X sin 9 A cos uj At) 



where the relative angular velocity Wb — \^-b ^ ~ ^d- 
In general, these relations represent anharmonically oscil- 
lating forces but we are taking R/d ^ 1 and they then 
reduce to an approximate form which is harmonic with fre- 
quency uj-B ~ f^d- The vertical force oscillates around its 
mean value with amplitude 



3{GA'h/d'^)iR/d) cost 



; sm ( 
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and the radial force oscillates with amplitude 
w (GMB/d^)sin6lB. 

3. Magnitudes of the neutron-star asymmetries 

In this section we give estimates for the magnitudes of 
neutron-star asymmetries arising in different ways. We first 
consider classical crystalline mountains and magnetically- 
confined accretion columns; both of these are found to be 
inadequate for the present purposes, however. We then 
turn to some different observationally-motivated possibil- 
ities which seem to be more promising. 

3.1. Isolated crystalline mountains 

Assuming that the basic nature of a mountain on the sur- 
face of a neutron star is the same as for mountains on plan- 
ets, the pressure at the base of the mountain needs to be 
less than the maximum shear stress of the surface mate- 
rial. This pressure is given by Pmnt = Pmntgns^mnt, where 
Pmnt is the average density of the material in the moun- 
tain, g^s is the surface gravity of the neutron star and /imnt 
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is the height of the mountain. The base of the mountain 
would be located at the outermost solid surface layer of 
the neutron star. The relevant density to take for this layer 
is rather uncertain; we will take it as being ~ 10^ g cm~^ 
and put Pmnt equal to that. The surface gravity is given 
by 5ns = GMp^/R\. For a neutron star of mass 1.4 Mq and 
radius 10km, q,is is 1.87 x lO^'^ c ms"^ . 

Following IStrohmaver et al] (jl99l[ ). the shear modulus 
of the neutron star surface material is taken to be 



0.1194 



M = 



n{Zef 



1 + 1.781 X (lGG/r)2 



(10) 



where n is the number density of ions, a is the inter- 
ionic distance, Z is the atomic number of the dominant 
ionic species, and F is the Coulomb coupling parameter 
(F > 10'^ for all practical purposes). The maximum shear 
strain in the surface of th e neutron star has been calcu - 
lated to be e ~ IQ-^-lQ-^ (|Smoluchowski fc WelcEflOTOh . 
although the re are suggestions that it might also be as 
high as 10~^ (jUshomirskv et all . l2000h . Taking typical val- 
ues Z = 26, n ~ lO^^cm-^ and a ~ 7 x 10"^" cm, the 
corresponding maximum shear stress is then S = iiO ~ 
10^^-10^^ dyncm"^. The maximum height of the mountain 
is obtained by setting Pmnt = S, which gives ^^n? i^ the 
range 0.01-lOcm . Taking the highest of these values, we 
then get the maximum possible mass of the mountain as 
being 

m„,nt - Pmnt(/iZ?)' 10^ g ^ lO-^^-lQ-^S 

which is too small to be relevant here. 



(11) 



3.2. Accretion columns 

For a neutron star with a strong magnetic field, accreting 
matter close to it can be diverted away from the equato- 
rial pl ane and form accret i on columns above the magn etic 
poles (jWooslev fc Wallace! . Il982t iHameurv et all I1983D . If 
the amount of matter in the columns is sufficiently large, 
this can provide another source of neutron-star asymme- 
try which could again be modelled in terms of "moun- 
tains" (probably two symmetric ones in this case). This has 
been invoked in connection w ith gravitational-wave emis- 
sion ([Melatos fc Pavn^ 120051 ) and we check here whether 
it could also be relevant in the present context. 

The column height can be determined from the condi- 
tion that the flow will start to spread out sideways when 
the pressure of the matter in the column Pac becomes large 
enough to bend the magnetic field lines outwards, typically 
when it is about a hundred times greater than the confining 
magnetic pressure, i.e. when 



Pac ~ 4 x lO^'^dyncm"^ (Ps/10^^ G) 



(12) 



where Bg is the s treng th of the surface dipole field 
(|Brown fc Bildstenl . Il998f ). For having hydrostatic equi- 
librium, this pressure should be the same as that else- 
where at the same level in the crust. In the density range 
lO^gcm"'^ < p < 10^°gcm~^ the pressure can be ex- 
pressed using the fit ting formula logP — 13.65 -I- 1.45 logp 
(|Bavm et al.l Il97lt) and the relation between the field 
strength and the density at the bottom of the column 
is then given by (pbot/10^ g cm'^) ~ 36{BJ10^'^ G)^-^^ . 
Using this, the scale height of the column, h^c is then 

Kc Pac/(pbot 9ns) ^ lO^cm (P./IO^^ G)" ''' . (13) 



Following IShapiro fc Teukolskvl (|l983h . the cross-sectional 
area of the column is estimated as 



10^° cm^ 



Ps 



1012 G 
R 



-4/7 



1.4M, 



1/7 



10^ cm 



9/7 



M 



2/7 



lO-SMo/yr 



,(14) 



where M is the accretion rate, which we normalise to a 
typical value for LMXBs. The mass of the column is then 
given by 



ma 



Pbot^xs^a 



Ps 



1012 G 
Rj, 



10/7 



1.4M, 



-6/7 



10^ cm 



9/7 



2/7 



lO-^Mo/yr 



(15) 



which is again very small, even for a field as high as lO^^ G. 
The continued existence of a roughly Keplerian accretion 
disc down to small radii, as required for our QPO pic- 
ture, needs the magnetic field to be suitably low and so 
magnetically-confined accretion columns do not seem to be 
relevant for the present purposes. 

3.3. Quadrupole moments inferred from limiting neutron star 
spin rates 

iBildstenI (|l998l ) noted that the spin frequencies for many 
accreting weakly-magnetised neutron stars were thought 
to lie in a rather narrow range around 300 Hz, which is 
a much lower frequency than that corresponding to cen- 
trifugal break-up. Since these objects are thought to have 
been accreting for long enough so as to gain sufficient angu- 
lar momentum to reach the break-up limit, it seemed that 
some mechanism was halting the spin-up. Despite the fact 
that the spr ead of spin frequencies is now thought to be less 
peaked (see [Manchester et all . 120051 ) , the issue of why spin- 
up should be halted before the break-up limit is reached 
remains a relevant one. Bildsten suggested that this might 
be caused by the accretion torque becoming balanced by a 
gravitational-wave torque resulting from an asymmetry of 
the neutron star. Taking the asymmetry to be represented 
by an Z = m = 2 perturbation, one finds that the mag- 
nitude of the misaligned quadrupole moment required for 
attaining this equilibrium at a frequency is given by 



)eq = 3.5 X 10^7 gcm2 ( JH± 

^ \1.4:MqJ VlO^cm 

1/2 



1/4 



/ Pa 



1/4 



.,_4t-)' (^V'"". (16) 



10-9Mo/yr/ V 300 Hz 



In terms of our simplified model of representing the asym- 
metry by means of point masses on the surface of an oth- 
erwise spherical neutron star, this corresponds to 



TOquad ~ 10 ^A/a 



1.4 M© 
M 



-3/4 



Pa 



106 ( 



-7/4 



1/2 



lO-SMo/yr 



300 Hz 



-5/2 



(17) 



Z. Stuchlik et al.: Gravitational excitation of QPOs (RN) 



4 



Various mechanisms have been suggested in the hterature 
for producing such a value (or even higher) , involving global 
deformations of the neu tron star rather than an isolated 
mountain on the sur face (lBildstenl . ll998HUshomirskv et alj . 
I2000t iHaskell et all . 120061, etc.). We will not enter into de- 
tails of this here but we note that values of the deformation 
large enough to be capable of explaining the limitation of 
neutron star spin-up, as suggested by Bildsten, may also 
be large enough for the present purposes: in the next sec- 
tion, we take the value for the point mass given by Eq. (fT7|) 
with the canonical parameter values, and check whether 
this same value ~ 10~^ A/a might also be large enough to 
be relevant for inducing the QPO behaviour. 



a plausible mechanism for inducing kHz QPO behaviour 
although many details remain to be worked out (in partic- 
ular concerning the response of the fluid medium and the 
production of the luminosity variations). The influence of 
the binary companion could possibly be effective in provid- 
ing the perturbations but the more likely possibility is that 
they might be produced by the neutron-star asymmetries. 
It is striking that the same magnitude for the misaligned 
quadrupole moment as advocated elsewhere for explaining 
limiting neutron star periods, seems also to give a plausible 
mechanism for inducing QPO behaviour. 



4. Discussion and conclusions 

Here, we discuss whether the effect of the neutron-star 
asymmetries or of the binary companion could be large 
enough to account for excitation of the QPO phenomenon. 
In our simplified picture (c.f. Eq. (H))), the excitation time 
for the amplitude a of resonant oscillations to grow to a 
particular value is given by 



(18) 



where /p is the amplitude of the perturbing force acting on 
unit mass, /o = GM^/R^ is the main gravitational force 
from the central object, Tk = 2^/^!^^ is the period of cir- 
cular Keplerian motion at the location being considered 
and the epicyclic frequency of the perturbation being ex- 
cited is a; = ar^K- The dimensionless radial and vertical 
"epicyclic functions" satisfy a < 1 everywhere. Note that 
since Tk 10^"^ s, the amplitude amplification in 1 second 
is ~ 10'^. The ratio a/R needs to grow to > 10^"^ in order to 
potentially explain the QPO behaviour and it would need 
to do that within ~ lO'^ s to account for the QPO phenom- 
ena seen in atoll sources. 

In the case of a binary companion, for the radial per- 
turbing force one has 



fo 



(MA (?: 



10" 



0.1 M» 



R 10^° cm 
106 cm d 



(19) 



where we have normalised to typical parameter values. This 
could produce a/R ^ 10^'^ at times t^x ^ lO'^s, but there 
is a problem for it producing resonances in the innermost 
parts of the disc because of having Wb ~ ^^k- The corre- 
sponding vertical force is smaller by at least four orders of 
magnitude and is clearly irrelevant. (Note that the verti- 
cal force is a tidal force whereas the radial one is a direct 
gravitational attraction; also, 6^ is probably rather close to 
7r/2.) 

For the neutron star asymmetries, we focus on the case 
of the misaligned quadrupole moments where /p//o can be 
~ 10~* in the inner parts of the disc for both radial and 
vertical oscillations (taking mquad ~ 10^^ Ma). For this, 
a/R could reach 10"'^ in < 10^ s, which encourages further 
investigation of this scenario. 

We conclude that at least one of the types of gravita- 
tional perturbation considered in this paper might provide 
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